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Twisted mass fermions allow light quarks to be explored but with the consequence that there are 
mass splittings, such as between the neutral and charged pion. Using a direct calculation of the 
connected neutral pion correlator and stochastic methods to evaluate the disconnected correla- 
tions, we determine the neutral pion mass. We explore the dependence on lattice spacing and 
quark mass in quenched QCD. For dynamical QCD, we determine the sign of the splitting which 
is linked, via chiral PT, to the nature of the phase transition at small quark mass. 



XXIII International Symposium on Lattice Field Theory 
25th-30th July 2005 

School of Mathematics, Trinity College, Dublin 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa . it / 



TM disconnected 



J. Pickavance 



1. TMQCD - Introduction 

TMQCD has Wilson quarks with the mass term modified by m —> m + ipLT^Y*,. In the contin- 
uum limit, a chiral rotation shows that this is equivalent to conventional QCD. At non-zero lattice 
spacing, however, the twisted mass formulation has several advantages, especially at maximal twist 
(where the effective mass term is imaginary). These advantages include (see ^ for a review) 

• improved behaviour of inverters for small pion mass, 



• order a 2 discretisation errors of flavour-conserving quantities, but with, however, the disad- 
vantage that there will be flavour breaking effects (eg. m(% + ) ^ ra(7T )) at non-zero lattice 
spacing 

Here we explore the n + , 71° mass difference both in quenched QCD to check scaling (with 
a) and in full QCD (Nf = 2) to explore the relationship via Chiral PT models to the nature of the 
phase transition at small quark masses. 

At maximal (90°) twist, the operators appropriate to create hadrons are shown below: 
hadron operator Disconnected piece 
k + uy 5 d 

71° ulu + dld Re Tr G u 

T]° ulu — dld ilm Tr G u 
a* uld 

c?q uY 5 u + dj 5 d ReTry 5 G„ 

fo uy 5 u - dy 5 d ilm Tr y 5 G„ 

where we have made use of the identity for propagators that G c \ = In several cases we need 

disconnected contributions: the notation is that TrG u (?) = X^M(jc,/)u(;c,f) an d the disconnected 
contribution is the correlation of G u (0)G u (t), with the vacuum piece of G u subtracted. This is 
illustrated by the right hand contribution below: 






2. Lattice evaluation of 7t° mass 

To evaluate the disconnected contributions to the hadronic correlation, we use a stochastic 
volume source with enough samples (24) that the stochastic sampling error is not dominant. We 
measure all qTq loops at all t (for zero momentum) and use both local and fuzzed sources. We 
use variance reduction [||] (gaussian random noise with subtraction of the lowest-order hopping 
parameter term) and we check the stochastic method by also evaluating the connected 7T° correlator 
which we compare with conventional (fixed source) evaluation. 

For quenched lattices at j3 = 6.0 and 5.85 with Wilson gauge action, K was chosen ^ to 
minimise parity mixing at /J. = 0. The quenched results presented here are also reported in ref [0]. 

As shown in fig. la, we find that the disconnected contribution is dominant in the n° correlator. 
This is advantageous since the large signal is more easily evaluated using stochastic methods. 
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Figure 1: (a) Ratio of pion correlators (with local operators) for quenched fermions; (b) Ratio of discon- 
nected to connected contribution to 7t° for dynamical fermions 



Table 1: Pion masses (in lattice units) 



N/ 


N g 


j3 


K 


a/» 


r /a 


71+ 


7T°(conn) 


7r°(tot) 


71° — 71+ 





400 


6.0 


.157409 


.0038 


5.368 


.1217(66) 


.218(9) 


.19(2) 


.07(2) 





400 


6.0 


.157409 


.0076 


5.368 


.1708(50) 


.246(5) 


•19(2) 


.02(2) 





400 


6.0 


.157409 


.0109 


5.368 


.2047(27) 


.267(4) 


.23(2) 


.02(2) 





100 


5.85 


.162379 


.0050 


4.067 


.1640(23) 


.319(14) 


.28(5) 


.12(5) 





100 


5.85 


.162379 


.0100 


4.067 


.2289(17) 


.351(8) 


.25(5) 


.02(5) 





100 


5.85 


.162379 


.0144 


4.067 


.2736(14) 


.380(6) 


.28(6) 


.00(6) 


2 


250 


.74 


.1590 


.0075 


3.87 


.1954(22) 


.37(4) 


.24(3) 


.04(3) 


2 


225 


3.9 


.1610 


.0075 


5.28 


.214(7) 


.35(2) 


.27(2) 


.06(2) 



We also plot the ratio of the full K correlator to the %+ correlator, and this is not consistent 
with 1.0. We quantify this flavour splitting by fitting both 71° and 7i + correlators separately (using 
both local and non-local sources) and we show our results in Table 1 . 

As well as evaluating the full tcq correlator, we can take advantage of the more accurate results 
available for the connected contribution to the 71° by interpreting [||] this contribution as from a 
(quenched) theory with fermion action with mass term m + iixy 5 . That theory has no disconnected 
contribution to TCo (since isospin is now a good symmetry) but the difference between it and the 
TMQCD value for the n+ mass can be used to explore lattice artifacts. These results are illustrated 
in fig. 2. Here we see that the 7T° - %+ mass difference decreases consistently with a 2 as a — ► 0. 
Moreover, we find that we can summarise the splitting by rQ(m 2 (7T°) — m 2 (7t + )) « c(a/ro) 2 , with 
the difference smaller including the disconnected contribution to the 71° (c « 10 but with large 
error) than connected alone(c « 23 and quite well determined). 

We also explore the n° mass with dynamical fermions [Q] - using Nf = 2 flavours of quarks 
(u and d). We fitted the 2x2 matrix of local and fuzzed correlators where both connected and 
disconnected contributions were here evaluated stochastically. We concentrate on results with ap- 



033/3 



TM disconnected 



J. Pickavance 




0.00 0.02 0.04 0.06 0.08 0.00 0.03 0.04 0.06 0.08 

(a/r ) 2 (a/r ) a 

Figure 2: Mass splittings between n° and K + in quenched studies. 

proximately maximal twist. Using the DBW2 gauge action (j3 = 0.74) (see also fig. lb), the pion 
masses are recorded in Table 1. The mass splitting m{n ,) —m(n + ) is consistent with zero but 
with a sign which corresponds (via leading order % PT [||]) to an Aoki phase at small quark mass. 
We also analyse results from the tree-level improved Symanzik gauge action (TLS j8 = 3.9) which 
are also presented in Table 1. The case nearest to maximal twist again shows a mass splitting 
m(7l°) — m(n + ) with a sign which corresponds with similar assumptions to an Aoki phase. 

3. Other mesons 

The flavour singlet pseudoscalar (T] for Nf = 2) involves a disconnected contribution even for 
conventional Wilson fermions. Here we evaluate the contribution for our twisted data sets. The 
full contribution is C — D where C is connected and D is disconnected. Positivity then requires 
C - D > or, equivalently, D/C < 1. 

We see, in fig. 3a that this requirement is violated in quenched QCD. Quenched QCD has 
ghost states and the flavour-singlet sector is not physical. For Nf = 2, however, the requirement 
D/C < 1 is well satisfied and it would be possible to extract the v\ mass from the rate at which 
it approaches 1, albeit with large errors. This is confirmation that the algorithms used to evaluate 
dynamical gauge configurations are consistent with positivity requirements. 

We also explore the spectrum of scalar mesons. In non-twisted QCD, the flavour non-singlet 
ao receives only connected contributions while the singlet fo receives both connected and discon- 
nected. In TMQCD, the ao states are again split with the a$ having only connected contributions 
whereas ajj has both. As shown in fig. 4, the disconnected contribution to a[j does help to make it 
closer to the charged ao, although there is a residual flavour violation (that we find decreases as the 
lattice spacing decreases). 

Again, because this time of the decay to r\ %, the quenched correlators are unphysical (they are 
negative). The Nf = 2 case, however, shows that positivity is restored, as it should be. 
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Figure 3: Flavour singlet pseudoscalar meson: ratio of disconnected to connected contribution to correla- 
tions. Here the total is C — D, so positivity implies D/C < 1. (a) quenched and (b) with Nf = 2. 



In detail, the scalar meson contributions will have to be corrected for the small admixture of 
pseudoscalar meson arising if the mixing angle is not exactly 90°. This small admixture becomes 
significant at larger t because the pseudoscalar mass is lighter than the scalar mass. 
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Figure 4: Scalar meson contributions (a) quenched and (b) with Nf = 2 



4. Summary 

• The disconnected contribution to 71° is big compared to the connected. 

• The 7T + 71° mass difference decreases with decreasing a as r^m 2 ^ ) — m 2 (7i + )) m c(a/ro) 2 
in a quenched study with c s» 10. 

• From a full-QCD study, one can extract the chiral lagrangian parameter relevant to 7T 7T + 
splitting. For the DBW2 gauge action, the splitting is small (consistent with zero within 



033/5 



TM disconnected 



J. Pickavance 



errors) with a sign which corresponds to an Aoki phase at small quark mass. For the TLS 
gauge action, the same sign is also seen. 

• The flavour singlet pseudoscalar meson signal (T]) is sensible and useful for Nf = 2 with 
dynamical TMQCD. 

• For scalar mesons, the need for disconnected diagrams and the possibility of the signal being 
polluted by mixing (if twist not 90°) with pseudoscalar suggests this will be a difficult area. 
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